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We explore the use of ‘optimized’ operators, designed to interpolate only a single meson eigenstate, 
in three-point correlation functions with a vector-current insertion. These operators are constructed 
as linear combinations in a large basis of meson interpolating fields using a variational analysis of 
matrices of two-point correlation functions. After performing such a determination at both zero 
and non-zero momentum, we compute three-point functions and are able to study radiative transi¬ 
tion matrix elements featuring excited state mesons. The required two- and three-point correlation 
functions are efficiently computed using the distillation framework in which there is a factorization 
between quark propagation and operator construction, allowing for a large number of meson oper¬ 
ators of definite momentum to be considered. We illustrate the method with a calculation using 
anisotopic lattices having three flavors of dynamical quark all tuned to the physical strange quark 
mass, considering form-factors and transitions of pseudoscalar and vector meson excitations. The 
dependence on photon virtuality for a number of form-factors and transitions is extracted and some 
discussion of excited-state phenomenology is presented. 

PACS numbers: f2.38.Gc, 13.40.Gp, 13.40.Hq, f4.40.Be 


I. INTRODUCTION 

The coupling of mesons and baryons to photons 
at leading order in a em is given by matrix ele¬ 
ments of the quark-field vector current (/|j 7 ‘|*) where 
j^{x) = i/>(x) 7 In the case that the initial and fi¬ 
nal state hadron are the same we express the matrix el¬ 
ement in terms of Lorentz-invariant form-factors, whose 
dependence on the virtuality of the photon, Q 2 , can be 
related to quark charge and current distributions within 
the hadron. The current can also induce a transition from 
one hadron eigenstate to another | /), in which case we 
speak of transition form-factors. For hadrons with non¬ 
zero spin, there are multiple possible amplitudes which 
can be labeled by the helicity of the hadrons, or we may 
expand the current in terms of multipoles to provide an¬ 
other convenient physically motivated basis. 

In the meson sector, these matrix elements appear in 
photo- and electro-production of mesons from nucleon 
and nuclear targets, where the coupling of the photon 
to a t-channel meson exchange is described by transi¬ 
tion form-factors. Measurements of these processes with 
unprecedented statistics will be made in the GlueX and 
CLAS12 detectors at the 12 GeV upgraded CEBAF [T|. 
In particular, photoproduction has been proposed as a 
means to produce large numbers of exotic J pc hybrids, 
those mesons which contain an excitation of the gluonic 
field as well as the usual quark-antiquark pair 0 , H- 
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Production of excited mesons at very forward angles 
using pion beams is also driven by the vector-current 
transition form-factor. In the Primakoff process the pion 
absorbs a nearly on-shell photon from a nucleon or nu¬ 
clear target, with an amplitude to transition to another 
meson species given by transition form-factors at Q 2 « 0. 
Recent such measurements of the couplings <22 —> 727 and 
7 T 2 -a 7 ry have been made at COMPASS j-J]. 

In charmonium and bottomonium, the relatively small 
total widths of the low-lying states means that radia¬ 
tive transition rates between them constitute significant 
branching fractions, and can be measured directly, as can 
rates of decay to a photon plus light-quark mesons where 
the heavy quark-antiquark pair annihilates 00 - 

In the baryon sector, the Q 2 dependence of transition 
form-factors can be measured quite directly in electro¬ 
production of excited nucleons off proton and neutron 
targets. The relative magnitudes of the various multi¬ 
pole or helicity amplitudes and the variation with pho¬ 
ton virtuality have been discussed as a means to study 
the internal quark-gluon structure of excited N* and A* 
states [7]. 

The properties of hadrons constructed from strongly 
interacting quarks and gluons should be calculable within 
the relevant gauge field theory, Quantum Chromodynam¬ 
ics (QCD). At the energy scale of hadrons, QCD does not 
have a small coupling constant and must be treated non- 
perturbatively. The tool we will use to achieve this is 
lattice QCD, in which the field theory is discretised on a 
finite grid of Euclidean space-time points, and where we 
can compute correlation functions as an average over a 
finite but large number of possible gauge-field configura¬ 
tions. 

The vector current matrix elements that we are in- 
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terested in, <( / 1 J M | i), appear in three-point correlation 
functions of generic form, (0\O f(tf) j^(t) O\(ti)\oy Here 
O], Oj r are operators constructed from quark and gluon 
fields, capable of interpolating mesons from the vacuum. 
In general, such hadron operators with definite quantum 
numbers interpolate not just one QCD eigenstate, but 
rather a linear superposition of all states with those quan¬ 
tum numbers. Each state propagates through Euclidean 
time with a factor e~ Et such that at large times only the 
state of lowest energy survives, suggesting that if we sep¬ 
arate the three operators in the correlation function by 
large time intervals we will obtain the transition between 
the lightest states with the quantum numbers of Oi, Of. 
However in many cases we actually wish to study states 
which are not the lightest with a given set of quantum 
numbers. If we use a generic operator to interpolate them 
from the vacuum their contribution will be as sublead¬ 
ing exponential dependences in the correlation function, 
which become dominated by lighter states as we prop¬ 
agate in Euclidean time. Determining the amplitude of 
subleading exponential contributions through fitting the 
time-dependence proves to be unreliable. 

Our solution to this problem is to form ‘optimal’ op¬ 
erators for interpolation of each state that we wish to 
study, that is operators which have a dominant ampli¬ 
tude to produce a particular state, and significantly re¬ 
duced amplitude to produce all other states, particularly 
those lighter than the state in question. In this case the 
three-point correlation functions are dominated by the 
desired initial and final states, even if they are not the 
lightest in a given quantum number sector. Optimal op¬ 
erators can be constructed as a linear superposition of 
operators in a large basis, where there is a particular lin¬ 
ear superposition for each state in the tower of excited 
states. The different superpositions are orthogonal in a 
suitable sense. We are able to find appropriate super¬ 
positions through a variational analysis of the matrix of 
two-point correlation functions, ( 0 |C , i (t)C ) ]( 0 )| 0 ), for a 
set of operators {Oi}. 

The technology we will utilize, constructing ‘opti¬ 
mized’ operators as linear superpositions of a large ba¬ 
sis of interpolators, is useful not only for extracting 
transitions featuring excited states, but also to improve 
ground-state signals by significantly reducing the un¬ 
wanted contributions of excited states to correlation func¬ 
tions. The presence of such contributions remains a prob¬ 
lem for calculations attempting precision extraction of 
ground-state matrix elements mm- Increase in statis¬ 
tical noise precludes separating operators by large time 
separations, so instead the use of ‘optimized’ operators 
seeks to deal with the problem directly by suppressing 
creation of the unwanted excited states. 

The use of optimized operators as a tool to extract 
excited-state radiative transitions was previously consid¬ 
ered in [15] , where the meson source operators, 0\ were 
optimized, while the meson sink operators, Of, were 
simple local fermion bilinears. Quark propagation from 


the meson sink proceeded using the ‘sequential source’ 
method, which proved to be a significant limitation on 
what could be achieved. In this study we will make use 
of the distillation framework for correlator construction 
m within which we will be able to use optimized oper¬ 
ators of definite momentum at both source and sink as 
well as to insert a vector current operator of definite mo¬ 
mentum. This is the first use of distillation to compute 
three-point functions. 

In order to demonstrate the technology we perform a 
calculation on dynamical lattices having three flavors of 
quark all tuned to approximately the physical strange 
quark mass - the spectrum of isovector mesons on these 
lattices was previously presented in [ia m- We con¬ 
sider pseudoscalar and vector mesons in the SU(3) flavor 
octet, extracting form-factors and transition form-factors 
for both ground and excited states in these channels. 
This should be expected to be a challenging undertak¬ 
ing - transitions between excited-state and ground-state 
vector and pseudoscalar mesons are ‘hindered’ and have 
relatively small magnitude for small photon virtualities 
- as such we will be extracting small signals from corre¬ 
lation functions built using ‘optimal’ excited state oper¬ 
ators. 

In section [TT] we introduce the decomposition of vector 
current matrix elements in terms of Lorentz covariant 
kinematic factors multiplying the unknown form-factors 
we seek to extract. We proceed to introduce the gauge 
configurations, operator basis, and background informa¬ 
tion relevant to two-point calculations in Sections |III| and 
IV Section [V] concerns three-point functions; we describe 


their calculation within the distillation framework and 
their relation to the matrix elements of interest while 
also demonstrating the efficacy of optimized operators. 


We present the results of our calculation in Section VI 


comparing with relevant previous calculations, and then 
conclude with a summary and outlook in Section [VII| Ap¬ 
pendices describing helicity operators, momentum con¬ 
servation in a finite-volume, the improvement of the vec¬ 
tor current and the flavor structure of the current follow. 


II. FORM-FACTORS AND TRANSITIONS 

The photon couples to the electric charges 
of u, d, s quarks via the vector current 
= +| u^u — | d'f^d — | sy^s, up to a factor of 
the magnitude of the electron charge, e. In general a 
transition induced by this current between a hadron, h , 
of spin-J and a hadron h' of spin- J' is described by the 
matrix element, 


where the spin-state of h is specified in terms of its he¬ 
licity, A, the projection of J along the direction of mo¬ 
mentum p. These matrix elements are simply related 
to the helicity amplitude for the transition 7 h —» h' by 
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including the initial-state photon’s polarization vector, 


M(^{X 1 ,q) hj(X,p) -» h'j,(X,p ')) 

= ^{X^,q)(h' J ,{X',p')\j> i \hj{X,p)), 

where q = p' — p and where the photon has a virtuality 

Q 2 = -q 2 = \p' -p| 2 - ( E h '(p ') - E h (p)) 2 . 

There are relations between these amplitudes which 
follow from the constraints of Lorentz invariance, cur¬ 
rent conservation and invariance under parity transfor¬ 
mations. These can be accounted for if we write a matrix- 
element decomposition in terms of a number of Lorentz 
invariant form-factors, Fi(Q 2 ), 


(h'j,(X! ,p')\f\hj(X,p)) 

= [hji (X',p'); hj(X,p)\ Fi(Q 2 ). (1) 


The Lorentz covariant ‘kinematic factors’, ivf, are con¬ 
structed from the meson four-momenta, p v , p’ v , and ini¬ 
tial and final state polarization tensors relevant to the 
spin of the mesons, e pCT . ..(A,p), e* a (X\p'). For any 
given pair of mesons h, h!. of definite spin and parity, 
there are only a limited number of possible construc¬ 
tions consistent with parity invariance and with the ad¬ 
ditional constraint of current conservation, we can write 
explicit decompositions in terms of a few independent 
form-factors. 

For example, a pseudoscalar particle like the pion has 
only a single form-factor appearing in its decomposition, 

{n + {p')\f\n + {pj) = {p+p'Y F n (Q 2 ). (2) 


At Q 2 = 0, the vector current measures the charge of the 
pion in units of e, so F„{ 0) = 1 exactly. 

In a transition between two different pseudoscalar par¬ 
ticles, there is again only one form-factor, but the kine¬ 
matic factor differs owing to the differing masses (to, in') 
of the pseudoscalar particles (n,n'), 


(n' + {p')\f\n + {p)) 


{p+p'Y+ 1j ^q^{p'-pY f^(Q 2 )- ( 3 ) 


The transition matrix-element between a vector parti¬ 
cle and a pseudoscalar can be expressed as 

{k + (p')\f\p + (X,p)) 

= & vP °p' v p p e g (A,p) m7r l mp F P ^(Q 2 ), (4) 

and for a vector meson stable under the strong 
interactions, the transition form-factor at Q 2 = 0 
can be related to the radiative decay width 
r(p+ 7T+7) = |F’ p7r (0)| 2 , where q is the 

momentum of the final-state photon in the rest-frame of 
the decaying p meson. 

The vector current matrix element for a stable vector 
hadron of mass to has three possible covariant structures 


having the right parity transformation properties once 
current conservation is demanded [T9] . 

{p + {X' ,p')\f\p + {X,p)) 

= - [{p+p'Y e*(A ',p')-e(X,p)\ G\{Q 2 ) 

+ [e M ( X,p)e*(X’,p')-p+e li *(X',p')e{X,p)-p '] G 2 {Q 2 ) 
- [{p + p'Y e*(X,p')‘P z{\p)-p' 2^\ G 3 {Q 2 ), (5) 

with a corresponding set of three independent dimension¬ 
less form-factors G\, G 2 , G 3 . A convenient basis having a 
clearer physical motivation is provided by the expansion 
of the vector current in terms of multipoles [20] . which 
in this case leads to a set of form-factors, 



G m = G 2 


Gq — Gi — G 2 + (l + J G 3, ( 6 ) 

which are proportional to the charge (Go), magnetic 
dipole {Mi), and quadrupole (G 2 ) multipoles respec- 
tivehQ At Q 2 = 0 they are related to the charge, mag¬ 
netic moment and quadrupole moment of the vector me¬ 
son: Gc{ 0) = 1, G m ( 0) = 2 to • p p , Gq{ 0) = to 2 • Q p . 

The other form-factors we considered above may also 
be identified with a particular multipolarity - in the 
p —> 7 T 7 transition case the single form-factor is of mag¬ 
netic dipole (Mi) type, while for the 7 r cases it is a charge 
form-factor (Go). For stable meson states, time-reversal 
invariance indicates that the form-factors are real func¬ 
tions of Q 2 . 

III. CALCULATION DETAILS 

In this first investigation of the extraction of excited 
state form-factors using distillation, we restrict ourselves 
to a single ensemble of gauge-field configurations, having 
three degenerate flavors of dynamical quarks tuned to 
approximately the physical strange quark mass. This 
set of anisotropic Clover lattices !22] (2B] has been used 
previously in studies of the meson spectrum mmm- 
[26] , meson decay constants [27] , baryon spectrum [28M3T] 
and meson-meson scattering [52H55] . For the calculations 
reported on in this paper, we used 535 configurations 
of lattice volume ( L/a s ) 3 x ( T/a t ) = 16 3 x 128, with 
a spatial grid spacing of a s ~ 0.12 fm and a temporal 
spacing roughly 3.5 times smaller. 

In this calculation we have an exact 51/(3) flavor sym¬ 
metry such that all the octet mesons ( 7 r, K, 77 ) are de¬ 
generate with a mass close to 700 MeV. Where results 
are expressed in dimensionful units, they are determined 


1 Note that the relationship between G 2 and Gm was presented 

with a typographic error in eh which is corrected here. 
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from the dimensionless quantities atE using the scale¬ 
setting procedure, 


E = 


atE 

a t mn 


m 


phys. 

n 


where a t ?no is the Q baryon mass calculated on this lat¬ 
tice and WQ llys ' is the experimental value [5]. 

Our use of a Clover-improved anisotropic quark action 
introduces an improvement term into the vector current 
which appears at tree-level. Discussion of the effect of 
improvement and renormalisation of the vector current 
will appear in Section [YE} 


IV. TWO-POINT FUNCTION ANALYSIS 

In order to determine radiative transition amplitudes 
between meson states within QCD we must first obtain 
the spectrum of states and find operators, constructed 
from quark and gluon fields, that reliably interpolate the 
states of interest from the vacuum. In general a color- 
singlet operator 0\ having definite J PC can produce all 
QCD eigenstates having those quantum numbers, 

o li°> = Ej n >< n i°li 0 >- 

We seek to determine optimized interpolators, fin, which 
when acting on the vacuum strongly interpolate only a 
single state with much reduced contributions from other 
states, 

fi n|0> = |n)(n|nt|0) + Y |m)(m|nj t |0) 

m/n 

= |n)(n|ftt| 0 ) + Y |m)e m - 

m/n 

In essence we seek a procedure by which we can minimize 
the e m (m 7 ^ n) relative to the strength with which our 
operator creates the n’th state, (n|fi„| 0 ). 

We will proceed by using a basis of interpolators, {Oi}, 
to construct two-point correlation functions of the form, 

where operators Oi are color-singlet constructions built 
from the basic quark and gluon fields of QCD, having the 
quantum numbers of the desired hadrons. Such correla¬ 
tion functions can be expressed as, 


A. Variational Analysis 

We propose that within any basis of operators there is a 
particular linear combination that is most suited to inter¬ 
polate the lightest state of the spectrum, another linear 
combination that optimally interpolates the first excited 
state, a third combination for the second excited state 
and so on. Thus optimized interpolators take the form 
fin = J2i ^”^ 1 ; where one can show that the best esti¬ 
mate for the weights w^ n \ in a variational sense, comes 
from solving the generalized eigenvalue problem |361l39j . 

C(t)vW=\ n (t)C(to)vW. (7) 

Here C(t) is the N x N matrix whose elements are the 
correlation functions Gjj (t) constructed from the basis 
of N operators, {Oi}, and is a generalized eigen¬ 
vector. The generalized eigenvalues, or principal corre¬ 
lators, A n (t) behave like e ~ E <'( t ~ t o ) a t large times, and 
can be used to determine the spectrum of energy eigen¬ 
states. The vectors are orthogonal on a metric, 
i>(m)tC'(to)u( n ) = <$ mitj where to is a reference timeslice. 
Examination of the orthogonality condition suggests that 
to should be chosen to be sufficiently large such that 
the correlation functions are dominated by the N lowest- 
lying states, with heavier states having decayed exponen¬ 
tially to a negligible level. Further considerations on the 
choice of t 0 are presented in [5HJ :35] ■ 

In practice we solve Eqn [7] independently on each 
timeslice, t, so that for each state, n, we obtain a time 
series of generalized eigenvectors to), which we ob¬ 

serve to be essentially time-independent with a suitably 
large choice of to- In practice we use the mean values 
(over the ensemble of gauge configurations) of the ele¬ 
ments chosen on a single timeslice to construct the 
optimized operators as 

nt = y^e-JW 2 Y v^ol (8) 

where the coefficients multiplying the sum are chosen to 
give the normalization (n|fl„| 0 ) = 2 E n , which will prove 
to be convenient when considering three-point functions. 

The procedure of solving the generalized eigenvalue 
problem in a basis of interpolating fields is carried out 
independently for each quantum number channel at each 
possible allowed momentum value. 


Cij{t) = Y 

n 


-^< 0 | 0 ? ( 0 )|n)<n|C»)( 0 )| 0 )e-^ i , 


where the spectrum of eigenstates is seen to control the 
Euclidean time dependence^ 
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We have introduced our particular choice of state normalization 
in a finite-volume here - we discuss this in Appendix Ib| 


B. Meson operator construction 

A straightforward approach to constructing a basis of 
operators capable of interpolating meson states is to use 
fermion bilinears containing some number of spatially di¬ 
rected gauge-covariant derivatives, that is operators of 
generic structure, 

O ~ • • • 1ft ip, 


(9) 






5 


where D = D-D. By expressing the vector-like deriva¬ 
tives and Gamma matrices in a circular basis, we can eas¬ 
ily construct operators of definite spin using the standard 
SO( 3) Clebsch-Gordon coefficients [T7] [TSj. Operators of 
this type with non-zero momentum can be constructed 
to have definite helicity, as described in [40]. 

In this calculation, QCD is discretised on a grid of 
points whose spatial structure has a cubic symmetry and 
as such we have not the complete continuous rotational 
symmetry of the continuum, but rather a reduced sym¬ 
metry: the symmetry of the cubic group at rest and the 
relevant little group in a moving frame. A consequence is 
that instead of having an infinite number of irreducible 
representations labeled by integer spin J (at rest) we only 
have access to the finite number of irreducible representa¬ 
tions, or irreps , of the cube, labelled Ai, T), T 2 , E, and 
A 2 . A corresponding argument applies in-flight, where 
the continuum helicity labeling is broken down to a finite 
number of irreps of the little group, see HD] for details. 

The operators of definite J (or helicity) constructed 
above can be projected into irreps of the relevant sym¬ 
metry group using a procedure called subduction , 

° 1 aI= E SifO JM (10) 

M=~J 

where A labels the cubic irrep and fi is the ‘row’ of the 
irrep (/j, = l...dim(A)). The subduction coefficients, 
Stf are tabulated in 18, 40;. These operators have been 
used extensively to study the excited spectrum of mesons 

ini omasa [3i]- 

C. Correlator construction through distillation 

Operators which interpolate hadrons from the vacuum 
have long been known to do so more effectively if the 
quark fields are suitably smeared over space H3 Eg. 
An extremely convenient method to do this is provided 
by distillation m, where the smearing operator is con¬ 
structed on time-slice t as an outer product of vectors in 
color and x— space, 

No 

□^w = E4 n) w4 n)t W’ ( n ) 

n—1 

where the Np vectors should be constructed to have 
strong overlap onto the low-energy quark modes most 
relevant to low-lying hadron states. A suitable choice for 
the vectors are the eigenvectors of the gauge-covariant 
three-dimensional Laplacian on a time-slice ordered by 
their eigenvalue. 

Smearing each quark field, a meson creation operator 
of fermion bilinear form with momentum p is 

<D\p) = ^tPxy e~ iM Tgz □jSfV’iS 

where time, color, and spin indices have been sup¬ 
pressed for brevity and where repeated position indices 


are summed. The object T can be non-local in x-space, 
and may for example feature gauge-covariant derivatives 
as discussed in the previous section. 

An advantage of the distillation framework is that it 
leads to a factorization of two-point functions into matri¬ 
ces describing quark propagation, called perambulators , 
and matrices describing operator construction. A generic 
connected two-point function using fermion bilinear con¬ 
structions, in which we explicitly show the smearing op¬ 
erators, can be decomposed as 

(o|'i/inr / n^(t) ^□ro^(o)|o) 

= ~Tnm( 0, t) $4p0) T pq(t> °) &qn (°) 

where the perambulators, T pq (t , 0) = ^ p ^(t) M t “ q 1 £^(0), 
can be obtained by inverting the Dirac matrix M 
on sources u (9) } ? = 1 .. .No a t timeslice 0, and where 
(t) encodes the operator construc¬ 
tion (here we include the momentum projection in T). 

D. Meson spectra & optimized operators 

We computed correlation matrices for irreps corre¬ 
sponding to J p = 0 _ , 1" at rest, and to magnitude of he¬ 
licity | A| =0,1 with nonzero momentum. We made use of 
np = [0,0,0], [0,0,1], [0,1,1], [1,1,1], [0,0,2], where the 
momentum is expressed in units of 2tt / L, p = (2tt/L) ftp. 
The quark flavor constructions were chosen to give ac¬ 
cess to the members of the SU(3)f octet - in this case 
the two-point function Wick contraction contains only a 
single connected diagram. 

In the rest frame we used operator constructions of 
the type in Eq. [9] with up to three derivatives while for 
nonzero momentum we used up to two-derivative con¬ 
structions. Further details about the construction of 
derivative operators and the basis used in this calcula¬ 
tion can be found in [T7j [18], 26, 40]. 

The resulting correlation matrices were analyzed using 
Eq.0 with each principal correlator, A n (£), being fit with 
the form, 

A„(t) = (1 - A n ) e+ A n e-^ 4 - 4 ”), (12) 

where the fit parameters are E n , E' n , and A n - the second 
exponential is present to absorb the effect of any states 
other than |n) remaining in the principal correlator. In 
practice, for suitably large values of to, we find that the 
energy scale of E' n is typically at or above the largest 
energy extracted, -E n =dim(C) ■ 

A typical example is presented in Figure [l] where the 
principal correlators for the lightest three states in the 
np = [0,0,0], A PC = Tj (A]~ + ) irreps, which contain 
1 (0 h ) mesons, are shown. We make use of a basis 

of 26 operators in the T l channel and 12 operators for 
A —*~ 

Another example is shown in Figure [2f a) for the 
case np = [0,0,1], A c = A]", which contains the he¬ 
licity zero component of mesons of ‘unnatural parity’ 
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FIG. 1. Left(Right) columns: Principal correlators, plotted as e E " ( ' t ~ t ° ) X n (t), for lightest three states in irrep .4j~ + (T 1 __ ). 
Two exponential fits shown, with the resulting mass spectrum of 0 ^ and 1 mesons shown in the central column. Radiative 
transitions and form-factors to be presented in this paper shown by the lines joining the states (solid lines are ‘charge’ transitions, 
dashed lines are ‘magnetic’ transitions). 


( J p = 0 _ , 1 + , 2 — ,...), where we present the effective 
mass of the principal correlator, m e g = ^ log 
with St = 3 at- As will be discussed in the next sub¬ 
section, these states appear to correspond (in order of 
increasing energy) to the ground-state 7r, the ai, the first- 
excited 7r', and the 7T2. 

We note that in both figures we observe that there is 
negligible curvature present for t > to, indicating that a 
single state is dominating the correlation function. This 
observation only holds for sufficiently large choices of to, 
and is one guide in the selection of a suitable to value. 
The time beyond which excited state contributions are 
negligible, which we can call t*, plays an important role 
in the construction of our three-point functions in terms 
of selecting the time separation of the source and sink 
meson operators. We desire our three-point functions to 
feature a time region in which the correlation function 
is dominated by the transition of interest, with contri¬ 
butions from other states having decayed away. In order 
to achieve this we should separate the source and sink 
projected operators (Eq. |8| by at least t* + tj. It re¬ 
mains possible that the vector current insertion may act 


to suppress or amplify the contribution of unwanted ex¬ 
cited states, but since we do not have this information 
in advance of the calculation, the above time separation 
serves as a reasonable estimation of the minimum. 

The source and sink cannot be separated arbitrarily 
far as the statistical noise on the entire three-point corre¬ 
lation function grows exponentially with increasing sep¬ 
aration. We can obtain an estimate of the maximum 
practical time separation for three-point functions by ex¬ 
amining the growth of noise in two-point function prin¬ 
cipal correlators. As an example in Figure [l] we see that 
the ground-state signal remains of high statistical qual¬ 
ity out to 25 timeslices (and beyond), while the first and 
second excited states begin to show significant fluctua¬ 
tions above t/a t ~ 15. Later we will find that while 
optimized ground-state three-point correlation functions 
are still statistically precise for time separations as high 
as 36 timeslices, excited-state correlation functions are 
not well determined for separations larger than around 
20 timeslices. 

In practice we solve the matrix problem, Eq. [TJ inde¬ 
pendently on each timeslice (as described in BSD- For 
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FIG. 2. (a) Effective masses of principal correlators for 

four lightest states in the irrep for momentum direction 
riff = [1, 0, 0] along with the energy determined from a two ex¬ 
ponential fit. (b) ‘Optimized’ operator correlation functions, 
(2E n ) _1 e Bnt (n n (t) ff£(0)), for the four states shown above. 
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FIG. 3. The rest-frame pion correlation function us¬ 
ing the (distillation smeared) i/rysi / operator (red) ver¬ 
sus using the ‘optimized’ operator, flo (blue). Plotted is 
2m. e m ^ <0|C*(£)C >t (0)|0> / |<0 |C|tt>| 2 . 
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sufficiently large to we find that the elements of the 
eigenvectors so obtained, v^ n \t), are essentially flat for 
t > to, and in practice we construct our projected op¬ 
erators, Eq. [8j using ensemble mean values taken from 
a single timeslice. In Figure [2](b) we show the corre¬ 
lation functions (0|fl n (t)fln(0)|0), observing that they 
behave in the manner we expect for optimized opera¬ 
tors. The corresponding off-diagonal correlation func¬ 
tions (0|Sl n (t)flJ l ,(0)|0) for n ^ W are statistically com¬ 


patible with zero for t > to. 

We conclude this section by demonstrating that an op¬ 
timized ground-state operator constructed as described 
above, does indeed significantly reduce the contribution 
of excited states to a correlation function, relaxing to 
the ground-state more quickly. This can be seen in Fig¬ 
ure^ where the optimized ground-state operator in A/ + 
is compared with the (distillation smeared) ip'Ysip opera¬ 
tor. 


E. Meson dispersion relations 


As previously mentioned, we independently compute 
the energy spectrum of states at each allowed value of 
p for each relevant irrep, but we expect to see the same 
mesons appearing at each momentum with an energy de¬ 
termined by their rest mass and the relevant dispersion 
relation. After identifying the mesons at each momen¬ 
tum (by their overlap with characteristic operators 140] ). 
we may examine their dispersion relation, E(\p |). This is 
presented in Figure [4] for mesons, 7r, p, p' and mesons 
cio, ai, bi,TTi (not used in this analysis) where they are all 
observed to be compatible with the relativistic dispersion 
relation, E 2 = m 2 + p 2 , or in temporal lattice units, 

(a t E) 2 = (a t m) 2 + (j^y) \np\ 2 , 

for a meson of mass m, with momentum 
citp = ^a s p = j/j^-[n x ,n y ,n z ], with the anisotropy 
taking the value £ = 3.44. Making use of optimized 
operators with momenta up to [0,0,2] allows us to 
sample many values of Q 2 in the form-factor extraction. 
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FIG. 4. Squared energies as a function of \riff\ 2 = (^) 2 |p| 2 for 
selected meson states. The points represent extracted energies 
for all irreps considered in the variational analysis (e.g. a 
vector meson like the p appears in three irreps with [1,1,0]: 
Ai, Bi, B?), while the lines show the relativistic dispersion 
relation for an anisotropy £ = 3.44. 
























V. THREE-POINT FUNCTIONS 


At 


t 


0 


We now turn to the three-point correlation functions 
used in this analysis which contain the vector-current ma¬ 
trix elements of interest. Their basic form is 


CV(Af,f) = <O|0 f (At) j p (t) 0\ (0)10), 

where the operators Oit are capable of interpolating me¬ 
son states of definite momentum from the vacuum - a 
suitable basis was discussed in the previous section. The 
relation between the correlation function and the desired 
matrix-element is exposed by inserting complete sets of 
eigenstates with the quantum numbers of 0\ and Of, and 
evolving all operators back to the origin of Euclidean 
time, 


C m (At,t) = 


ni,rif 


1 i 

2E nt 2 E n . 


g E nf (Ai t) g E n A 


X <0 1 O f (0) | nf ) (ttf |(0) | tii) (tii 1 0\ (0) 1 0). 


The summation runs over all states, but clearly if the sep¬ 
arations between the operators are large, At t 0, 
only the lightest states in the i and / channels will con¬ 
tribute and we can extract the vector-current matrix el¬ 
ement between them. However, at modest time sepa¬ 
rations there will remain subleading exponential contri¬ 
butions from excited states, and these ‘polluting’ terms 
can be a source of systematic error in the extraction of 
ground state matrix elements [THUiS]. Reducing ex¬ 
cited state pollution by simply separating operators by 
longer Euclidean times is not always practical, due to the 
increase in statistical noise with increasing separation. 

One of our major aims here is to extract excited-state 
matrix elements, which we may access using the opti¬ 
mized operators described in the previous section. Using 
the optimized operator for an excited state should lead to 
a three-point correlation function whose leading behavior 
at large times is given not by the ground state, but rather 
by the relevant excited state. If we are interested in the 
ground-state there is also an advantage to using the ap¬ 
propriate optimized operator in that it will have much re¬ 
duced overlap onto low-lying excitations (relative to any 
single operator in the original basis, for example V’TsVh 
c.f. Figure [3]), leading to a corresponding reduction in 
the excited state pollution in the three-point correlator. 

Three-point correlation functions using optimized op¬ 
erators take the form, 

Un f fini (At,t) = (o|n nf (At)j'*(i)nt 1 (o)|o) 

= e ^M At - 4) e- E 'i* (tif 1^(0) |rii) + ... 

(13) 


where the leading time-dependence is that of the states 
(|rii), |rif)), selected by the choice of optimized operators. 
The ellipsis represents the residual contributions of other 
states, which should be significantly suppressed when us¬ 
ing optimized operators - we will explore the degree to 
which this is manifested in explicit calculation. 
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FIG. 5. Graphical depiction of a connected three-point cor¬ 
relator. The gray blobs represent the source and sink oper¬ 
ators carrying the quantum numbers of the initial and final 
state mesons. The blue line corresponds to a perambulator 
(T pq (At,0) = while the red line rep¬ 

resents a generalized perambulator as described in the text, 
Qp q (At,t,0) = f (p) t(At)M^ t 1 t r t M t 7 o 1 ^ (,) (0), which carries 
the momentum and quantum numbers of the current inser¬ 
tion. 


A. Correlator construction & distillation 


Three-point correlation functions featuring a current 
insertion require a slight extension of the distillation 
framework presented in Section |IV C| since the quark 
fields in the current should be those which appear in the 
action, which are not smeared. Exposing the smearing 
operators, the general form required is, 

■ ^□rav>(0)|0), 


which, considering for now only the completely connected 
Wick contraction, can be decomposed as 


(o|^nr / nv>(At) • ip(t) ■ '^/’□r*n-?/j(o)|o^ 

= -W0, A t) <f> f mp (At) [^(AiJM^rM^^O)] $‘„(0) 
= —T„ m (0, At) $^ P (A t) g% q (At,t,0) $® n (0), 


where the object in square brackets, Q pq (At, t, 0), is a 
generalized perambulator. It can be obtained through in¬ 
version from sources {£^} 9 =i...at d at timeslice 0 to ob¬ 
tain 0), inversion from sources {£^} P =i...jv D at 

timeslice At to obtain M t “ ^j£00(A t) which can be re¬ 
lated to £( p )t(Ausing q 5 hermiticity, and con¬ 
traction with the operator insertion, T, at each value of t 
between 0 and At. In this calculation we do not average 
over multiple time-sources separated by the same value 
of At, although this could be done to increase statistics. 

In our application, exposing the Dirac spin, 
x —space and color indices, the current insertion is 
^s^y-a.b = dz,yS a ,b- In the case of an improved 

vector current, which we consider later, we require also 
r £x?;a.& = a ik Sx,y8 a ,b, yielding another set of gen¬ 

eralized perambulators. 

Within this construction we are able to project each 
operator into definite momentum, and as such we only 
compute correlation functions in which the momentum 
is conserved, Pf = p, + q. Some discussion of momentum 
conservation in a finite volume appears in Appendix [B| 
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B. Correlation functions using optimized operators 


Turning first to the case of three-point functions with 
pion-like operators at the source and sink, we plot in 
Figure [6] the form-factor (as defined in Eqn. [2]) extracted 
from the three-point function, 




where O represents either (in red) or the opti¬ 

mized operator f l n (in blue). The sink operator, located 
at At = 28 at ~ 0.9 fm, is in the A c = A£ irrep of mo¬ 
mentum np f = [1, 0,0], while the source operator, located 
at t = 0, is at rest in the A PC = A^ + irrep. We clearly 
observe that the optimized operators give rise to a signal 
which is flat over a number of timeslices away from the 
source and sink, corresponding to the contribution of just 
the ground-state pion, while the simpler operators 

over this time range always retain a non-negligible pollu¬ 
tion from excited states. Such behavior is expected from 
our two-point function analysis: for example, at rest we 

0.73, so the distillation smeared op- 


find 


(O|'075'0|i'i=l) 

<0|'i/'75V’|n=0> 


erator [V’Ts'i/'] , acting on the vacuum, creates both the 
ground and first excited state with comparable strength. 

Our principal motivation for using optimized operators 
is to get access to transitions involving excited states. 
In Figure [7] we show matrix elements extracted from 
three-point correlation functions computed using either 
the ground-state 7r or first-excited state 7 r' optimized op¬ 
erator at the source (t* = 0, p, = [-1, 0,-1]) and either 
the ground-state p or first-excited state p' operator at 
the sink (tf = 20 a tl Pf = [1,0,-1])- We observe that 
there are clear statistically significant signals for excited- 
state transitions when using the appropriate optimized 
operators. 

In general, even for optimized operators, there may still 
be some residual contamination coming from states that 


2 



lie beyond the reach of our variational basis, and indeed 
curvature away from flat behavior as we approach the 
source or sink timeslice is observed in Figures [6] and [7J 

In order to make maximal use of the time-series data, 
in particular in those regions where there remains some 
unwanted excited-state contribution, we opt to perform 



0 2 4 6 8 10 12 14 16 18 20 


FIG. 6. Form-factor from vector current three-point function 
with pion operators at source ( t = 0, n = [0,0,0]) and sink 
(At = 28 at, n = [1,0,0]). Red points correspond to using 
the ‘unoptimized’ bilinear i/jysV'; the ‘optimized’ operator is 
shown in blue. 


FIG. 7. Form-factors extracted from optimized three-point 
correlation functions with a 7r or 7r' operator with p = [-1, 0, -1] 
at t = 0 and a p or p' operator with p = [1,0, -1] at t = 20 at. 
The source-sink separation in physical units is roughly 0.7 fm. 
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FIG. 8. Upper panel shows the pion ground-state form-factor extracted from correlation functions with source-sink separations 
of At = (12,16, 28) at using the procedure of fitting with Eqn. 14 


rif. 


= [1,0,0], rtp f = [2,0,0], 


The lower panel illustrates this using the example of 


a correlated fit over a time range with the form, 

F(Q 2 -1) = F(Q 2 ) + f f e~ SEf ( At ~^ + ft e~ 5Ei * (14) 

where //, SEf , /), SEi and F(Q 2 ) are real ht parameters. 
We make further use only of the constant term, which 
corresponds to the desired form-factor. Fitting the data 
to this form also exposes the energy scale of the pollu¬ 
tion terms, SEf and SEi. Generically, when present, we 
find that these energies lie at or above the scale of the 
largest energies we reliably extract in our two-point func¬ 
tion variational analysis. In cases where there is a clear 
extended plateau region, we may exclude the exponential 
terms and perform a fit to a constant value. 

The dependence upon source-sink separation, At, for 
the ground-state pion form-factor can be seen in Figure[8] 
The lower panel shows the illustrative case ftp, = [1, 0, 0], 
np f = [2,0,0], where we observe that there is only a vis¬ 
ible plateau region for At = 28 at, while for the shorter 
separations, At = 12 at, 16 at, we make use of a fit using 
The resulting values of F n (Q 2 ) are observed 


14 


Eqn. 

to be compatible - other time separations were also ex¬ 
plored with similar results. The upper panel shows that 
this procedure is generally applicable and leads to form- 
factors from each time-separation that are in agreement 
across a range of Q 2 . 


In practice, while we extract a very large number of 
form-factor determinations at many Q 2 -values, we choose 
to make use of only those where application of Eq. [M] to 
F(Q 2 -,t) shows modest excited-state contributions. Any 
cases where a clear trend toward a constant value is not 
visible are discarded. 


C. Extracting multiple form-factors 

Eqn. |T] presents the general form of the decomposi¬ 
tion of a vector-current matrix element into independent 
form-factors, Fi(Q 2 ), and the corresponding kinematical 
factors, Ki, which depend upon momenta and helicities. 
Moving to a more complete notation including momen¬ 
tum and helicity labels, our three-point correlation func¬ 
tions may be written, 

<°l n n ( ,ft,A f (At)^(t)nt iiftiAi (0)|0) 

= e- E ''F At ~ t ' , e~ E ''i t '^2K?(ni,p { ,\[-,n i ,pi,\ i ) Fi{Q 2 ) + ..., 
i 

where as previously the ellipsis represents possible pol¬ 
lution from states other than (|ni), |tif)), which will have 
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residual time-dependence, but which as shown in the pre¬ 
vious section are suppressed when using optimized oper¬ 
ators. 

Any one correlator provides, in general, an under¬ 
determined linear system for the multiple form-factors 
we wish to extract. By using different combinations 
of initial and final state helicities and momenta all at 
the same Q 2 , we can build a linear system which is 
constrained or over-constrained, from which we can de¬ 
termine the set {Fi(Q 2 )}. Removing the known Eu¬ 
clidean time-dependence, e~ E "te~ E "i t , from the 
correlation functions, we should be left with objects 
which are time-independent up to pollution from other 
states, which should be modest for optimized opera¬ 
tors. For fixed state choices, n;, nf, using an indexing 
a = (pi, A;; n',pt, Af), we can write the linear system, 


r„(At,t) 

= e E 'S <0|fi nfiftiAf (A t)fpt) <, p ,, Ai (0)|0> 

= ]Ta pa) ^(Q 2 ) + ..., (15) 


which is of the form F = IK • F with F a vector over a, IK 
a rectangular matrix with indices a, i, and F a vector of 
form-factors, indexed by i. This may be converted into a 


system featuring a square matrix: IK"*" F = 


IKt IK 


F, which 


can be inverted using SVD. In the case where only a 
single form-factor contributes this procedure can still be 
followed as a way to average over rotationally equivalent 
momentum combinations. 

In practice we solve this system independently for each 
value of t between 0 and At if our optimized operators 
were perfect, leading to no pollution from other states, 
we would obtain the same form-factors on each timeslice. 
In fact we obtain Fi(Q 2 ;t), where the time-dependence 
is fitted as described in the previous section to account 
for the presence of pollution from other states. 

This approach was previously used in m and ng in 
the extraction of charmonium form-factors. 


D. Cubic symmetry 

A consequence of discretizing QCD on a hypercubic 
grid is that the theory does not possess the full three 
dimensional rotational symmetry of the continuum. In¬ 
stead, we are restricted to a subset of rotations that 
leave the cube invariant. This smaller symmetry group 
has only a finite number of irreducible representations 
into which the infinite set of continuum representations 
labelled by integer spin, J, must be subduced. A sim¬ 
ple example is J — 2, where the five equivalent ‘rows’ 
(M = —2... 2) get distributed into a three-dimensional 
irrep called T 2 and a two-dimensional irrep called E. Be¬ 
cause there are only a finite number of these irreps, they 
must accommodate multiple values of J , such that T 2 
also contains parts of J = 3,4.... For systems with non¬ 
zero momentum, the symmetry group is called the ‘little 


group’, and the corresponding subduction is from helic- 
ity, A. Tables of the spin/helicity content of cubic irreps 
can be found in (401. 


To correctly reflect the symmetry of our theory then, 
we should label our correlation functions according to ir¬ 
reducible representations of the cubic symmetry. In prac¬ 
tice this is what we do by computing using the subduced 
operators introduced in Section IV B Using these oper¬ 
ators, the three-point functions take the form, 


(o|q a ;;«(a <*^(0)10), (i6) 

where the indices A, /i label the cubic group irrep and 
the ‘row’ (1 • • • dim(A)) of the irrep. 

Considering only the cubic symmetry of the lattice, 
and not any underlying continuum-like symmetry, we 
would not expect there to be any relationship between 
different irreps. Furthermore matrix element decomposi¬ 
tions should be defined in terms of the irreps of the cube, 
not in terms of hadrons of definite spin. For example 
a correlation function with a T 2 operator at the source 
should take values that need not be related to one with 
an E operator at the source. 

However, were there really to be no relation, we could 
hardly claim to be approximating QCD in a realistic 
manner. In practical calculations it should be the case 
that through a combination of sufficiently fine lattice 
spacing, reduction of discretization artifacts through im¬ 
provement of the action [44j , and interpolation of hadrons 
using operators smoothed over many lattice sites JE5], 
that the continuum symmetry is manifested to a good 
approximation with only small deviations. For example 
we might expect to see a relation between X 2 and E corre¬ 
lation functions corresponding to them originating from 
the same J = 2 meson. In previous two-point function 
calculations we have observed that the rotational symme¬ 
try of the continuum theory is clearly visible in relations 
amongst the irreps both for eigenstate masses and the 
values of matrix elements (0|O^|n) TF, ITS] . 

Since we expect to see a comparable restoration of the 
rotational symmetry in this calculation, we do not at¬ 
tempt to build decompositions according to the symme¬ 
tries of the cube, rather making use of the continuum-like 
hclicity decompositions presented earlier, subduced into 
irreducible representations of the cube. 

A slight additional complication in this analysis arises 
from our use of anisotropic gauge configurations in 
which the space and time directions are discretized 
with different spacings. Spatially directed currents will 
need to be renormalized separately from temporal cur¬ 
rents and the discretization effects along the two direc¬ 
tions are expected to be different - in explicit calcu¬ 
lation we will not mix spatially directed currents with 
their temporal counterparts. Had we used isotropic 
lattices the temporal component of the vector cur¬ 
rent would be related to the spatial components, how¬ 
ever here we will keep them separate with the tem¬ 
poral component of the current subducing differently 
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FIG. 9. (a) The energy spectrum for the little-group irreps, A c = A 1 , B 1 , B 2 , A 2 with momentum rip = [1,1,0], which 

are observed to contain states with degeneracy pattern consistent with the lowest p(l~~) meson and the lowest bi(l " 1 ). 
(b) The transition form-factors extracted from the correlation functions (0|flp A ^ (28, p') j M (£, q) fit (0,p )|0) with npi = [0,1,1], 
rip = [0,-1,1] for the irreps A c = A ~[, Bj”, B 2 , which we observe to have consistent values, differing only in the amount of 
excited state pollution (the current is projected into the En irrep of momentum direction rip = [0, 2,0]). The black points in the 
upper left are the result of solving the linear system described in Section [V C[ including all equivalent rotations of the source 
and sink momenta. 


from the spatial components. For spatial component^] 


the subduced current is 


A.-y , /^-y 

> J= 1,A 


j x where 


j x = e(q, A) • j, whereas temporal components subduce as 

■A ryifl-y .p — fl 

0 


J 


In order to relate the irrep-based correlation functions 
that we compute, Eq. [l 6 j to the helicity-based decom¬ 
positions presented in Eq. |T] we define subduced matrix 
elements, which for the spatial current case take the form, 


<n f ,p f , A f , p { \j A r^-r |n i; pi, Ai, m) 

E qAf,fAf 

L D J 7 = 1 ,A 7 

Ai 5 A-y , Ai 




(17) 


A Y ' ^ 7 ) * Ki Pi)') ). 


A similar expression exists for the temporal portion of 
the current, here it will subduce into a “scalar” one di¬ 
mensional irrep. 

The use of the Lorentz-covariant decomposition in this 
expression implies relationships between different irreps 
that we must establish are present in the computed cor¬ 
relation functions for this approach to be considered rea¬ 
sonable. 


3 Here we have made the choice to treat the vector current as a 
creation operator as opposed to an annihilation operator. The 
alternate definition would induce changes of phases throughout 
the calculation (e.g. subduction coefficients and momentum pro¬ 
jection). 


In Figure ia) we show an example of the 
extracted spectrum across little-group irreps, 
A c = A^ ,B^ ,B 2 ,A 2 for np = [ 1 , 1 , 0 ], where the 
distribution of states matches the expected subduc¬ 
tion patterns for a pair of meson states, a lighter p 
( J PC = 1 ) state and a heavier b\ ( J PC = l +_ ) state. 

Forming the optimized operator for the p state in each of 
the A]”, B^ , B 2 irreps, we can compute the three-point 
function, 

(0|fi p (A q) 0^(0,p)|0), 


for ftp' = [ 0 , 1 , 1 ], ftp = [ 0 , - 1 , 1 ], and np = [ 0 , 2 , 0 ], across 
the little-group irreps A^,B^,B 2 at the sink. The 
source operator is the optimized operator for the ground- 
state pion in the A 2 irrep. The three different sink irrep 
choices correspond to the subduced versions of the three 
helicity projections of a vector meson. For At = 28at, 
the resulting form-factor is plotted in Figure [9^b) , where 
we observe that while the amount of excited state pollu¬ 
tion differs slightly in each irrep, the form-factor values 
are consistent, indicating that we are observing compo¬ 
nents of the same 1 meson in the three irreps. 

In general we find that the relationships between ir¬ 
reps implied by using Lorentz-covariant decompositions 


of matrix elements of hadrons of definite spin as in Eq. 17 


are present in our correlation functions and that the cor¬ 
responding linear systems of the type expressed in Eq. [T5] 
can be solved satisfactorily. 

The cubic nature of the boundary of the lattice has 
additional implications which impact the properties of 
unstable hadron states - we defer a discussion of this 
point to Section [VTT] 
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E. Renormalization & improvement of the vector 
current 


The local vector current is not conserved with a 

Clover discretized fermion action, and should be renor¬ 
malized multiplicatively by a factor Zy. We determine 
this factor non-perturbatively by computing the charge 
form-factor of the 7r + or p + meson at Q 2 = 0, where the 
continuum value corresponds to the charge of the meson 
in units of e, F( 0) = 1, so that 


F c ° nt .( 0 ) 1 

ww ~ ef 


(18) 


On an anisotropic lattice, where we have treated space 
and time differently in the action, there can be one 
Zy for the spatial vector current, ip'y z '4 , i and another 
for the temporal vector current, We extract 

the zero momentum transfer form-factor from correla¬ 
tion functions with identical source and sink momen¬ 
tum, (0|fl w (At,p) j^(t,q = 0) nj. (0,p)|0). In Figure [To 


we show our extracted values of Zy observing no signifi¬ 
cant dependence on the momentum p. We do however ob¬ 
serve some dependence upon whether we extract from the 
pion form-factor or from the p form-factor. The depen¬ 
dence upon the state, a discretization effect, is expected 
as we have not used a conserved current. In addition, 
some of the discrepancy can be attributed to imperfect 
tuning of the anisotropy parameters in the fermion and 
gauge action j22j , where the measured fermion anisotropy 
is 3.44 compared to the target value of 3.5. 

Using the n extraction, which is statistically most pre¬ 
cise, performing a correlated average over momenta, we 
obtain, 


= 0.846(6), Z$r = 0.961(7), (19) 

for the spatial and temporal renormalization factors re¬ 
spectively. All subsequent presentations of form-factor 


0 0.2 0.4 0.6 0.8 1.0 | p | 2 / GeV 2 



FIG. 10. Vector current renormalization factor extracted from 
Q 2 = 0 form-factors of the pion (circles) and the p (squares). 
Spatial current in blue, green, temporal current in red, orange. 


values in this paper have been multiplicatively renormal¬ 
ized by the ensemble mean value of the relevant choice 
of these two factor^ 

The anisotropic discretization also introduces a tree 
level O(a) improvement term not present in isotropic 
calculations which amounts to replacing the Euclidean 
current = '07 At ' i / ; with 


J4 = il’lii’ + \ y (! - 0 a sdj (ipa 4 jip) 

jk = + I (1 - £) a t d 4 (Vi(T4fcV’) 1 (20) 

where £ = a s /at = 3.44 is the anisotropy, and v s = 1.3 is 
a parameter appearing in the anisotropic fermion action 
[22l [46] . Details of the derivation of the improved current 
are deferred to Appendix |C| 

In Figure |TT] we plot our determination of the pion 
form-factor at a range of Q 2 values using both the unim¬ 
proved and the improved current. We note that the ad¬ 
dition of the improvement term, which over the Q 2 range 
considered provides only a small shift, does bring the spa¬ 
tial and temporal current extractions into better agree¬ 
ment. Since we observe the effect of the improvement 
term to be small, and our main aim is to explore the use 
of ‘optimized’ meson operators, in the remainder of this 
report, with the exception of the ground-state pion form- 
factor, we will make use of only the spatial component of 
the unimproved vector current. 
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FIG. 11. The pion ground- stat e form-factor for unimproved 
( , 07 m ' 0) an d improved (Eqn. 20) currents. 


4 the statistical uncertainty on Zy can be considered an overall 
systematic error on the normalization of form-factors. 
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VI. EXTRACTED FORM-FACTORS & 
TRANSITIONS 


In this section we present form-factors and tran¬ 
sitions for the lightest few isovector pseudoscalar 
and vector mesons. We make use of the current 
j v = +'^u r ] v u — ^drf v d — |s 7 ^s, such that the form- 
factors are in units of e, the magnitude of the elec¬ 
tron charge. This calculation is performed with three 
flavors of dynamical quark all having the same mass, 
tuned approximately to the physical strange quark mass. 
We extract vector current matrix elements between 
(/, I z ) = (1, +1) members of SU(3)f octets. Discon¬ 
nected diagrams do not contribute to the amplitudes con¬ 
sidered in this analysis as demonstrated in Appendix [D] 
where the flavor structure of the current is explored fur¬ 
ther. 


A. Form-factors 

1. 7r form-factor 

The pion form-factor appears in the matrix element 
decomposition, {^ + {p')\j^\^ + {p)) = {p + p'Y F V (Q 2 ), 
which we will extract from three-point Euclidean corre¬ 
lation functions computed using optimized ground-state 
pion operators of definite momentum at the source (at 
t = 0) and the sink (at At = 28 a t ). As discussed pre¬ 
viously, we will present F{Q 2 \t), where the leading Eu¬ 
clidean time-dependence of the correlation function has 
been removed, with any remaining time-dependence sig¬ 
naling the presence of excited state contributions to the 
correlation function. By utilizing many values of p and 
p' we can determine the form-factor at a range of Q 2 val¬ 
ues. We plot F^(Q 2 \ t) for a subset of these Q 2 values in 
Figure [l2j where for each Q 2 we overlay a fit according 
to the form in Eq. |T~T| 


0 0.2 0.4 0.6 0.8 t/fm F n (Q 2 ) 



FIG. 12. Typical F n (Q 2 ',t) extracted from optimized three 
point functions (points) with fit descriptions using Eqn. 
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(curves). Note that the data points have a high degree of 
timeslice correlation which is accounted for in the fitting. 


In Figure 13 we plot the resulting Q 2 dependence, 
shown via both dimensionless a 2 Q 2 and scale-set using 
the D-baryon mass prescription presented in Section m 
A large number of kinematic points are sampled by con¬ 
sidering all combinations of momentum such that n ?. < 4, 
rip, < 4 and < 4. The extracted points, for the im¬ 
proved current discussed in Section YE appear to lie 
on a single curve, with only small residual scatter which 
can originate from fitting-range systematics and modest 
discretization effects. 

Describing the Q 2 dependence may offer some phe¬ 
nomenological insight, albeit in this calculation at an 
unphysically heavy quark mass. A commonly used ap¬ 
proach to describe vector-current form-factors of hadrons 
is to argue that the photon is behaving like the lightest 
vector meson which can couple to the hadrons in ques¬ 
tion, which in this case would be the p. This “vector 
meson dominance” (VMD) describes the Q 2 dependence 
by EVmd(Q 2 ) = i + Q 2 / m 2 ■ Using the p mass determined 

on these lattices, m p = 1020(1) MeV, we have the dashed 
curve shown in Figure |T3j which is seen to describe the 
lattice data reasonably well only for small photon vir¬ 
tualities. One possible explanation of this effect is that 
as we move out to larger Q 2 , considering only the near¬ 
est time-like pole, the p, and neglecting all excitations, 
becomes a progressively poorer approximation. 

The distribution of charge within the pion can 
be characterized by the charge radius, defined via 
the slope of the form-factor at zero virtuality, 


>^6^F(Q 2 )| 


Q 2 =0‘ 


We may obtain this quan¬ 
tity from the discrete Q 2 data presented in Figure [l3] 
by parameterizing the Q 2 -dependence for small virtual¬ 
ities. Considering gaussian (F^(Q 2 ) = F(0)e~® A 6 / 3 " 


;) forms to de- 


and pole (F„(Q 2 ) = F( 0) 1+Q2/m2 , 
scribe <3 2 <0.3GeV 2 , we obtair0 a charge radius 
{r 2 )l /2 = 0.47(6) fm, where the error includes the vari¬ 
ation over fit-form. As we might expect, in a cal¬ 
culation where three flavors of quarks all have ap¬ 
proximately the strange quark mass, we obtain a pion 
charge radius somewhat smaller than the physical pion 
{r 2 )\ /2 = 0.67(1) fm (5, 37], and also smaller than the 
physical kaon (r 2 )^- 2 = 0.58(4) fm [48] . 


2. p form-factors 

The three form-factors required to describe the vector- 
current response of a vector hadron may be defined as 


5 If -F(O) is allowed to float in fits, a value statistically compatible 
with 1 is obtained, as it must since the pion form-factor at zero 
Q 2 was used to set Zy. The fit y 2 values obtained are fairly 
large due to the scatter in the statistically precise data, which is 
likely due to small discretization effects which are not described 
by these smooth fit-forms. 
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FIG. 13. Pion ground-state form-factor, F W (Q~), using the improved current, Eqn. 20 Vector meson dominance using the p 


meson mass on this lattice shown by the dashed curve. Fits to the small-Q 2 points using gaussian and single-pole forms shown 
by the gray curves. 


0.2 0.4 0.6 0.8 1.0 Q 2 / GeV 2 



FIG. 14. Ground-state p meson multipole form-factors. Points have the same color and shape labeling presented in Figure [T3| 
Fits to the Q 2 dependence, described in the text, are shown as gray curves. 


in Eq. |6j which makes use of a multipole basis. The de¬ 
composition presented in Eq. [5] defines the linear system 
which we may solve, as described in Section [V C for the 
form-factors. We plot the charge, Ge{Q 2 ), magnetic, 
Gm{Q 2 ), and quadrupole, Gq(Q 2 ) form-factors in Fig¬ 


ure [14] Examination of Eqs [5j [6] indicates that only the 
charge form-factor has a non-zero kinematic factor when 
Q 2 = 0, and as such only it is determined there, while all 
three form-factors are sampled for positive non-zero Q 2 . 
The smallest form-factor, Gq , shows the largest scatter, 
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which likely originates from modest discretization effects 
and timeslice fitting-range fluctuations. 

Fitting the Q 2 dependence of the charge form-factor 
with various fornisj^J over various Q 2 ranges we obtain 
Gc(0) = 0.94(1) and ( r 2 ) X J 2 = 0.55(5) fm where the er¬ 
rors include a systematic variation over different fit forms. 
The deviation of the charge from 1 was discussed previ¬ 
ously in Section |VE| 

In order to determine the magnetic and quadrupole 
moments from Gm( 0) and Gq( 0) it is necessary to 
parameterize the Q 2 dependence of the form-factors 
and extrapolate back to Q 2 = 0. Utilizing a range 
of possible forms, we obtain Gm(0) = 2.17(10) and 
Gq(0) = —0.54(10), accounting for the variation over fit- 
forms, which is much larger than the statistical uncer¬ 
tainty, in the errors. More precise determinations of these 
quantities could be obtained if twisted boundary condi¬ 
tions were used to sample the form-factors at smaller Q 2 
(see for example [49]). 

Within a simple picture of the p as a qq bound-state, 
the presence of a quadrupole moment would indicate a 
required admixture of D- wave into the dominantly S- 
wave wavefunction. Previous estimates of the /9-meson 
magnetic moment in versions of QCD with heavier than 
physical quarks come from chiral effective theory [50] 
where Gm{ 0) ~ 2.2 for large pion masses, and quenched 
lattice QCD using either an energy shift in a magnetic 
field f-TT] where Gm(0) = 2.13(6), or extrapolation to 
zero Q 2 from a single spacelike virtuality [52] where 
Gm( 0) = 2.05(4), at comparable unphysical pion masses. 
A dynamical calculation, Ref. [53] . which appeared while 
this manuscript was in the final stages of production, 
found, at a comparable pion mass, Gm(0) = 2.23(2) 
and Gq(0) = —0.362(20), using a model extrapolation 
to Q 2 = 0 from a single non-zero Q 2 point. 

3. 7r' form-factor 

The examples presented in the previous two subsec¬ 
tions were the lightest states with the relevant quantum 
numbers. As such it was not strictly necessary to use 
optimized operators - any suitable meson interpolators 
used in the three-point functions will, in the limit of large 
time separations, give access to the matrix elements. We 
will now move to the case of an excited state, the first 
excitation of the pion, which we access using optimized 
operators to eliminate the contribution of the ground- 
state pion. 

As described in Section[Vj the signals for excited states 
are typically noisier than those for the ground state, and 
as such we separate the source and sink operators by a 


6 G(Ole-^ 2 / 16 ^ 2 , G(o) e -e 2 ( 1 +“ < 3 2 )/ 16 ' 32 , G(0)/(1 + Q 2 /m 2 ), 
G(0)/(1 + Q 2 /m 2 + 7(Q 2 /m 2 ) 2 ), G(0) 


smaller time, in this case At = 16a*. The decomposition 
for this matrix element is of the same form as the pion 
described previously, Eq. [2] We plot the extracted form- 
factor, F^{Q 2 ]t), as a function of the current insertion 
timeslice in Figure [15] 

The Q 2 dependence of the form-factor, iv ( Q 2 ), is pre¬ 
sented in Figure[l6| While the extracted values at Q 2 = 0 
are not statistically precise, they are certainly consistent 
with unity. The charge radius can be extracted from the 
slope at Q 2 = 0 which we determine by parameterizing 
the data for Q 2 < 0.3 GeV 2 , yielding (r 2 )]/, 2 = 0.74(6) fm 
where the error includes variation over parameterization 
form. As we might expect for a state which likely can be 
characterized as a radial excitation, this is significantly 
larger than the 0.47(6) fm found for the ground-state pion 
at this quark mass. 

Ref. [53] , computing at a very similar pion mass found 
0.517(4) fm for the ground-state pion charge radius, and 
0.59(3) fm for the first excitation of the pion. Their ap¬ 
proach determines a single point on the form-factor curve 
at Q 2 ~ 0.16 GeV 2 which is used to determine the slope 
at Q 2 = 0 assuming monopole dependence on Q 2 . 

B. Radiative Transitions 

1. tt' try transition 

In a transition between different pseudoscalar mesons, 
the decomposition of the current in terms of a form-factor 
F^niQ 2 ) is as in Eqn.[3j and the form-factor must vanish 
at Q 2 = 0. The transition form-factor is extracted from 


0 0.1 0.2 0.3 0.4 0.5 t/fm J r -‘\Q 2 ) 



FIG. 15. Current insertion time dependence for the form- 
factor of the first excitation of the pion, F n /(Q 2 ;t), shown 
for a range of source and sink momenta. The high degree of 
timeslice-timeslice data correlation is manifested in the error 
on the fit which is not significantly reduced relative to the 
error on the individual data points. 


7 Gaussian (F w /(0)e Q 2 / 16 ^ 2 ) and one-pole (F^/fO)/!! + 
Q 2 /m 2 )) forms were used. 
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FIG. 16. First-excited pion form-factor, F n i (Q 2 )- Points have 
the same color and shape labeling presented in Figure |13| 
Fits to low Q 2 dependence used to constrain charge-radius, 
as described in text, shown as gray bands. 


three-point functions with At = 20 at, fitting the time- 
dependence as previously to account for any residual un¬ 
wanted excited state contribution. We plot the extracted 
form-factor in Figure 17 that we are now able to ex¬ 
plore the timelike Q 2 region, where previously all points 
were spacelike, follows from the differing masses of the 
hadrons at source and sink, a simple example being the 
case where p' = p , so that Q 2 = — ( E'(p) — E(p)) 2 < 0. 
In order to be able to trivially relate our Euclidean am¬ 
plitudes to Minkowski amplitudes, we must restrict our¬ 
selves to the region where the current is not timelike 
enough to produce on-shell hadrons. In this calculation 
where the 7T7 t threshold is above the p mass, this limits 
us to Q 2 > —m 2 ~ — 1 GeV 2 . In order to explore further 
into the timelike region, a somewhat more sophisticated 
approach must be followed (oT; 7 m] • 
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FIG. 17. Transition form-factor between first-excited and 
ground-state pions, F —> n 7 . Points have the same color 
and shape labeling presented in Figure [13] with the excited 
pion having momentum pf. 


2. p —► 7 T 7 transition 


A suitable decomposition for a vector to pseudoscalar 
transition in terms of a dimensionless form-factor is given 
in Eqn.[4j Using optimized operators for the ground-state 
p and ground-state 7 r we computed correlation functions 
with At = 28 at for a large range of source and sink mo¬ 
menta - the resulting determi nati on of the form-factor, 
Fpir(Q 2 ) is presented in Figure 
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The value of the form-factor at Q 2 = 0, known as the 
photocoupling, is of particular interest since it controls 
the rate of the physically allowed radiative tran sition pro¬ 
cess, p ± —» 7r ± 7. As can be seen in Figure [l8j we do not 
determine this quantity directly, but we may estimate it 
using interpolation between our space-like and time-like 
points. Using a range of fit forms over several Q 2 ranges 
(plotted in grey) we estimate F pn ( 0) = 0.494(8), where 
the error includes variation over fit-forms. 

The Lorentz invariant matrix element for the decay 
p + —► 7r + 7 can be obtained by contracting the matrix 
element in Eqn. [4] with a final state polarization vector, 
= e*(A 7 ,g)(7r + (p")|j M |/° + (A,p)), and for a vector 
stable under the strong interaction, we may obtain the 
decay width from 


I> + ^ ^ = 32^S2>^l a . 

where we have summed over the final state photon polar¬ 
izations and averaged over the initial state polarization 
of the p. Using the decomposition above, and restoring 
the factors of e, we obtain the result relating the width 
to the photocoupling, 


r(p+ —> 7r + 7) = a- 


IsT 


3 (rrip+m.ir) 2 


F P ,m 2 , 


where a — e 2 /47r. 

The calculation performed here uses three degenerate 
quark flavors tuned to approximate the physical strange 
quark mass and as such our photocoupling determination 
cannot be directly compared with experiment. For ori¬ 
entation we show in Figure [T8[ the experimental values, 
F pn { 0) = 0.33(2) and Fk*k( 0) = 0.57(3) extracted from 
the corresponding decay rates obtained via the Primakoff 
effect for pions and kaons incident on nuclear targets [5Sh 

ES|. 

The Q 2 dependence of this meson transition form- 
factor plays a role in models of deuteron electromagnetic 
structure, where a virtual photon probe may couple to 
the bound nucleons or to the meson currents proposed to 
supply the binding m- 


3. p' —> 777 transition 

The first-excited p state may also undergo a transition 
to the ground-state pion, with the form of the decompo- 










18 


0 0.2 0.4 0.6 0.8 1.0 Q 2 / GeV 2 



FIG. 18. Ground-state p to ground-state n transition form-factor. Curves in gray show fits used to interpolate between spacelike 
and timelike regions to determine the photocoupling, F pn (0). Experimental decay widths converted to photocouplings shown 
for orientation. 


sition of the matrix element being the same as in the pre¬ 
vious section. In Section |IV| we presented the spectrum 
of excited vector mesons, finding that the first-excited 
state, TO p / = 1882(11) MeV, is close to being degener¬ 
ate with the second-excited state m p " = 1992(6) MeV. 
Our use of optimized operators corresponding to orthog¬ 
onal combinations of basis operators allows us to reliably 
study the two excitations independently. 

We extract the form-factor using optimized operators 
in correlation functions with time-separation, At = 20 at, 
with the results presented in Figure [19} To determine the 
photocoupling, F p ^(0) = 0.050(4), we perform fits to the 
data over various Q 2 ranges using several fit-forms, and 
the quoted uncertainty includes this variation. 

The photocoupling for this transition is observed to 
be an order of magnitude smaller than that of p —> 7ry 
extracted in Section [Vi B 2| Within simple models treat¬ 
ing mesons as qq bound-states with non-relativistic wave- 
functions, such a suppression is expected - the net effect 
of the current is to slightly shift in momentum-space the 
wavefunction of the pion, and since the p' is likely de¬ 
scribed as a radial excitation, the resulting wavefunc¬ 
tion overlap is much reduced relative to that for the 
ground-state p. This is described as a ‘hindered’ mag¬ 
netic dipole transition. A relevant experimental example 
of a hindered transition lies in the charmonium sector 
- the relative rates of ip(2S) -7 r] c j and J/ip -7 ?y c 7 , 

r rS^ZyyWj/l ^ 1 ~ °’ 1 sh ° W the ex P ected hierarchy 
of hindered versus non-hindered |5]. 


p" —> Try transition 

An extraction analogous to that presented in the pre¬ 
vious subsection can be performed for the second-excited 
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0.10 

0.08 

^ 0.04 
0.02 


-0.02 -0.01 0 00 1 0.02 0 03 0 04 a 2 Q 2 

FIG. 19. First-excited p transition to ground-state 7r form- 
factor, Fp'^iQ 2 ). Points have the same color and shape la¬ 
beling presented in Figure [TS] Gray curves show fits used to 
interpolate to the photocoupling. 



p state, leading to the form-factor shown in Figure 20 
Interpolating to Q 2 = 0 using a range of forms yields 
Fp" v {Q 2 ) = —0.016(3), which is smaller still than the 
p' —> 7T7 photocoupling. The sign is somewhat arbitrary 
and would only have definite meaning were we to com¬ 
pare to other transitions involving the p". 


Within a simple qq bound-state model we might ex¬ 
pect the p" state to be dominated by a 3 D i configuration 
(and indeed the operator overlaps presented in [3] seem 
to suggest this), which would have a ‘hindered’ structure 
in a transition to the ground-state S-w&ve pseudoscalar 
owing to the need for the current to provide a D -wave 
angular dependence, which appears only as a relativistic 
correction to the leading behavior. 
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FIG. 20. Second-excited p transition to ground-state 7 r form- 
factor, F p n 7r (Q 2 ). Points have the same color and shape la¬ 
beling presented in Figure [18] Gray curves show fits used to 
interpolate to the photocoupling. 

5. n 1 —> p 7 transition 

The first-excited pion may undergo a transition to the 
ground-state p. The results, extracted from At = 20 a t 
correlation functions, are presented in Figure [2TJ along 
with a number of parameterizations used to interpolate 
a photocoupling of F 7 r / p (0) = 0.18(2). Again we observe 
a significant suppression relative to the p —> 777 case in 
line with this being a hindered transition. 
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FIG. 21. First-excited 7r transition to ground-state p, 
F 7T i p (Q 2 ). Gray curves show fits used to interpolate to the 
photocoupling. 


6. p' -7- 7/7 transition 

This transition, which occurs between excited states, 
is not expected to be hindered in the case that the p' and 
it' are identified predominantly as the first radial excita¬ 
tions of the p and n respectively. As such we might expect 
a somewhat larger photocoupling than in previous sub¬ 
sections. We extracted the form-factor from optimized 
operator correlation functions wit h At = 20 at obtaining 
the results presented in Figure [ 22 ] Fits to the Q 2 depen¬ 
dence with a range of forms lead to an estimate of the 
photocoupling, F p i n i(Q 2 ) = 0.7(2), which, although not 


determined with high precision, is of comparable size to 
the p —> ttj coupling. 


VII. SUMMARY 

A desire to extract current matrix elements between 
excited hadrons motivated our exploration of optimized 
operators which are capable of interpolating only a single 
hadron eigenstate from the vacuum rather than a super¬ 
position of all possible eigenstates. 

We have demonstrated through explicit calculation, 
the utility of these optimized operators when used at the 
source and sink of three-point correlation functions also 
featuring a vector current insertion. Matrix elements fea¬ 
turing the two lightest pseudoscalar mesons and the three 
lightest vector mesons were explored, with successful ex¬ 
traction of excited-state transition matrix elements at a 
range of momentum transfers. In the case that the op¬ 
timized operators correspond to ground-states, their use 
reduces the degree of unwanted excited-state contribu¬ 
tion to the correlation function and allows for the source 
and sink to be separated by a smaller time interval with 
a corresponding reduction in statistical noise. 

Optimized operators, as we have implemented them, 
are linear superpositions of a large basis of meson opera¬ 
tors, and it follows that we need a technique that allows 
efficient evaluation of potentially complicated operator 
constructions. Distillation has previously been shown to 
meet these needs in the case of two-point correlation func¬ 
tions, and in this study we have demonstrated its efficacy 
in the case of three-point functions, the first time it has 
been used in this manner. Its use allows the problem to 
be factorized into pieces corresponding to operator con¬ 
structions of definite momentum at source and sink and 
independent pieces corresponding to quark propagation. 
This separation of operators from propagators avoids the 
problem encountered in sequential-source techniques of 
requiring the operators to be selected before the Dirac 
matrix inversions take place. 

This first study was restricted to transitions involving 
pseudoscalar and vector mesons, but we may easily ex¬ 
tend it to other meson quantum numbers using the flex¬ 
ible operator basis presented in htuxsw One impor¬ 
tant application is the determination of transition pho¬ 
tocouplings relevant to meson photoproduction experi¬ 
ments like GlueX and CLAS12, where the process is mod¬ 
eled as proceeding through f-channel meson exchange. 
Part of the role of these experiments is to search for ex¬ 
otic and non-exotic hybrid mesons, those states which 
have an essential gluonic contribution in their wavefunc- 
tion - it has been claimed previously in models [2 [59] 
that transition photocouplings of hybrid mesons may be 
large, motivating photoproduction searches. The tech¬ 
nology we have explored for mesons may also be applied 
to the baryon sector in which the dependence of the tran¬ 
sition form-factors on photon virtuality can be measured 
in electroproduction experiments on proton and neutron 
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FIG. 22. First-excited 7r transition to first-excited p, F n i p i(Q 2 ). 


Gray curves show fits used to interpolate to the photocoupling. 


targets [TJ. 

Throughout this analysis we have proceeded under the 
assumption that the hadrons we consider are stable eigen¬ 
states of QCD. In fact, for the quark mass used in the 
calculation only the ground-state 7r and p are below all 
relevant kinematic thresholds, while the various excita¬ 
tions can in principle decay. As the light quark masses 
are reduced toward their physical value, even the ground 
state p ceases to be an eigenstate and rather appears as 
a resonance in tttt scattering. In recent years there has 
been significant progress determining meson resonance 
properties using lattice QCD, making use of the discrete 
spectrum of states in the finite volume defined by the pe¬ 
riodic lattice boundary, following the formalism initially 
presented in [60'. 

To date, however, there is no calculation exploring the 
coupling of a resonance to external currents. A rigor¬ 
ous calculation at physical kinematics, where the p is 
a resonance, seeking the coupling p —> 7ry would in fact 
need to determine the P -wave partial-wave amplitude for 
n'y —> tttt as a function of the invariant mass, m By 
analytically continuing the amplitude to complex values 
of and extrapolating to the p-resonance pole, the 
coupling could be extracted as the residue of the ampli¬ 
tude. Only very recently m has the formalism relating 
matrix elements extracted in a hnite-volume to the phys¬ 
ical amplitude been presented. 

The techniques laid out in this paper, which allow the 
extraction of matrix-elements for each state in a tower 
of discrete eigenstates, will be required in any attempt 
to determine resonance couplings to external currents. 
The extension of the operator basis to include, as well as 
single-meson-like operators, also meson-meson construc¬ 
tions has already been explored in two-point functions 
[55H55] . and the corresponding three-point function cal¬ 
culations can now be attempted. 
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Appendix A: Rotations and helicity operators 

In this analysis we have made use of sets of kinemat¬ 
ically equivalent helicity matrix elements. The opera¬ 
tors appearing at the source and sink of our three point 
functions are subduced, with the subduction coefficients 
implying certain choices of rotation conventions. One 
method to consistently relate kinematically equivalent 
subduced matrix elements to a canonical frame involves 
embedding additional phases into the correlation func¬ 
tions appearing in our linear system, Equation |15| Fol¬ 
lowing m a helicity state is defined by 

|p; J, A) = Rp L z (p )| J, A) (Al) 
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where | J, A) is a state at rest with spin J and J z com¬ 
ponent A, L z (p) is a boost along the 2 -axis with mo¬ 
mentum magnitude p , and Rp is a rotation that rotates 
from the 2 -axis to direction p. A set of conventions to 
implement these rotations, which separates them into a 
rotation from the 2 -axis to a reference momentum direc¬ 
tion, followed by a rotation which leaves the cubic lattice 
invariant is presented in (40] . where the corresponding 
subduction into irreducible representations of the little 
group is also shown. 

For an arbitrary rotation, R, these helicity states trans¬ 
form as 


R |p; J, A) = |Rfi J, A), (A2) 

where the helicity is left invariant and where rotations 
about the direction of the momentum, p, introduce a 
helicity-dependent phase. 

A canonical state , where the spin state is specified us¬ 
ing the projection along the 2 -axis is defined as 

|p; j , m) = Rp L z (p) Rp 1 | J, m) 

= L(p) | J,m), 

where L(p) is a Lorentz boost along direction p with mo¬ 
mentum magnitude, p. These states transform under ro¬ 
tations as 


helicity matrix elements in different fra mes. To achieve 
this we must determine the phases in Eq. 


A2 


0 i&(R,p,J, A) 


corresponding to the rotation conventions laid down in 
[401 used in the construction of our helicity operators. 

This can be achieved by evaluating matrix elements of the 
form (p'; J', A'| R \ p; J, A). It is convenient to exchange 
the helicity states for canonical states using Equation |A3| 

- the matrix element can then be written as 

(p'-,J',X'\R\p;J,X) 

= E D m'\* (ftp') E; m'\R\ p; j, m) D^l(Rp) 

m' ,m 

= Y, D ^y( R p') D ^ R ) D ^l( R p) 6 3 {p'- r p) 6j,j. 

m' ,m 


Using the composition of rotations we have 

Dy\{ R p' lRR p) = E D %X'( R P') D ™U R ) D ™l(Rp), 

m' ,m 

and we conclude that the phase in Eq. |A2| is given by 

e rHR,P, J +) = D w (R-lRRp) . (A4) 

A generic matrix-element of the vector current between 
states of definite momentum, spin and helicity, is 

(p 7 ; J', A'| | p; J, A), 


R\p;J,m) = J2 D m’m( R ) | R p;J,m'), 

m' 


and for an arbitrary rotation R this matrix element trans¬ 
forms as 


and are related to helicity states via 

|p;^,a) = E^mi^p) |p; j » to )> ( A3 ) 

m 


<P , ;^ , ,A'| |p; J,A) 

= (p'- 1 J',\'\R- 1 Rj> 1 R- 1 R\p-,J,\) 


where Rp is the rotation for direction p appearing in 

K<|.[AT] 

The application of this formalism that is required for 
the work reported on in this paper is to consistently relate 


since the current rotates as a four-vector, 
R R^ 1 = [ R ~*‘] l ^j v - Inserting complete sets of 
states, and using the fact, expressed in Eq. |A2| that 
rotations do not change the helicity of states, we obtain 


(p"; J', A'| f\ p; J, A) = [R~X Jd 3 q Jd 3 k (p 7 ; J', A'| R | g; J 7 , A')<<?; J 7 , X'\f \ *; J, A)(fc; J, A| R \ p;; J, A). 


(AS) 


Using our newly derived representation of the phase, Eq. Ad we can write this as 


(p'-J\X'\f | p; J, A) = [R~X Dx'\*(R~Rp’ RR &) D iR ( R Rl RR p)( R p'-,J',X'\r\Rp-,J,X) 


In practice we use this expression to relate matrix elements determined for various allowed lattice momenta, p, p ', 
to matrix elements for some reference momenta p re f = I?p, p r 7 ef = Rp'. 
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Appendix B: Momentum conservation in a finite-volume 


We define meson eigenstates which in infinite volume 
have normalization 

<n(fc)|n'(p)) = <5nn'(27 t) 3 2E n (k) S^ 3 \k-p), 


momenta of free particles is quantized, k = -^n^, 
where Hk = (n x ,n y ,n z ) and the normalization becomes 
(n(fc)|n'(p)) = 6 nn >2E n (k)L 3 and completeness 

is expressed as 


such that the completeness relation takes the form 



d 3 k 1 
(2tt) 3 2 E n (k) 


n(k))(n(k) 


In a periodic cubic volume, L x L x L, the allowed 


1 = -pr ^ —^^|n(fc))(n(fc)|. (Bl) 

L3 rtr 2E ^y 

Two-point correlation functions in which the source 
and sink operators are projected into definite momentum 
have a spectral representation which can be obtained by 
inserting Eq. |B1[ 


C(t) = (o\o { (p { ,t) Ol(pi,0)\0) 

= (°l o f (ar,i) 0)1°) 


= ^ZZ Jnrd L 


n n k 


2E n (k) 


ptk 


{L%%)e~ Ent <0|O f (x = 0,0)|n(/c))<n(fc)|O; (y = 0,0)|0) 


= ^f Eat <0|(0,0)|n(pi))<n(pi)|O/(0,0)|0), 


where we note an explicit factor of the lattice volume, L 3 . 

Three-point correlation functions projected into definite source, sink and current momentum have a spectra repre¬ 
sentation, 


C{t) = (0|O f (pf, tf) j(q, t) e>/(pi,ti)|0) 

1 1 


= L 3 S. 


'pt,Pi+q 


E 


2E n . 2E n 


e Bnf(t/ t] e B "i ( * u) (0|e> f (0,0)|n f (p f ))(n f (pf)|j(0,0)|n i (p i ))(n i (p i )|C> I t (0,0)|0>, 


which again features an explicit factor of the lattice volume, L 3 . This volume factor, common to two-point and 
three-point functions may conventionally be absorbed into the meson creation/annihilation matrix elements. 


Appendix C: Vector Current Improvement on 
Anisotropic Clover Lattices, 


The form of the anisotropic Clover fermion action that 
we use [23] 05] can be obtained from the desired Eu¬ 
clidean action using the field transformations 


(l T 2 atm T 2 Ut 7 4 T 2 ^)i ) 

T = V>(l + \Cl m a t m + 0(74^4 + if \ s a s 'j j ' < 5j) : 


leading to 




(1 + i(fl m + n m )a t m)m 

+ (l + t 0,tTTl) 74^4 + O'tE 4^4 

-f- ^1 —2 (Hm “l - flm)^#^^ T fls a s m) 7 j -t- Us a s l)j /),,■ 
“f~ 2 (fliUt + Hs^s)^i 4 T ,;4 T (fls^s) ^ ^ i.j 

i>j 


where a term proportional to {l?i, D 4 } was eliminated by 
choosing f2 t = —fit and fl s = —fls. Making the choices, 
fl m + fl m = 1 , fI s = — g^s, fit = — y, discretizing the 
derivatives, smearing in the spatial directions and insert¬ 
ing tadpole factors gives the action presented in [22] |25] . 
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Reference |22j describes the non-perturbative tuning of 
the parameters. 

The same transformation can be applied to the desired 
vector current to obtain the classically O(a) 

improved vector current. The transformation gives 

Jti = (l + \a t m) 

- - < 5 m4 ^(^4 - 

- \v s a s (dj{$cT^ip) - - ~dj)i/j ), 

and use of the classical equations of motion for the quark 
fields allows for the elimination of the gauge-covariant 
derivatives acting on quark fields to give 

J4 = (1 + \(m + m 0 )a t ) ^74 V’+ 27 (1 - £) a s dj 

jk = (l + + rn 0 £)at) ijj'yk'tp + 3(1 - 0 ntdi^a^ki’), 

where £ = a s /a t is the anisotropy. 

Since we non-perturbatively determine the renormal¬ 
ization of the vector current at any given quark mass by 
computing meson form-factors at zero momentum trans¬ 
fer, we choose to absorb the mass-dependent factor into 
the renormalization factor and consider the improved 
vector currents, 

J 4 = Z'y (^74+ 37K 1 - 0 a s dj (^(74^)) 

j k = Zy^ikip + j(l - 0 a t d4('ipcnki’)y (Cl) 

As expected we observe that the improvement terms van¬ 
ish at classical level in the case of an isotropic action, 

£ = i- 


Appendix D: Wick contractions for three quark 
flavors 

In this paper we work with a version of QCD in which 
the up, down and strange quarks all have the same mass, 
leading to an exact SU(3)p symmetry. We compute 
radiative transition matrix elements between states in 
octets, 8 p, of SU(3)f- In particular we consider the 
(I,I Z ) = ( 1 ,± 1 ) elements, and interpolate mesons from 
the vacuum with operators having flavor structure dTu. 

Formally integrating out fermions from the path- 
integral we replace field pairs q x q y with the appropriate 
quark propagator, Q xy , the inverse of the Dirac opera¬ 
tor. Since the Dirac operator depends upon the quark 
mass, there can be differing propagators for each flavor 
of quark, i.e. uu —► U, dd —»• D, ss —► S. In the case 
of exact isospin symmetry, m u = m d and U = D and if 
SU(3)f is exact, U = D = S. 

For transitions between 1=1 mesons, inserting cur¬ 
rents u'y f jii, d'y^d or sj u s, the schematic Wick contrac¬ 
tions shown in Figures [231(a) and (b) appear, where those 
shown in Figure |23[ b) are referred to as disconnected 
contributions. Were we to consider 1 = 0 mesons, we 



FIG. 23. Possible quark line contractions required for me¬ 
son three-point functions with a current insertion. Blue lines 
represent quark propagation which could be described by a 
perambulator and red lines represent quark propagation with 
a current insertion which could be described by a generalized 
perambulator. Green lines indicate quark propagation which 
does not begin or end on an operator which can be distillation 
smeared, and as such cannot be described by a perambulator. 


could also receive contributions from the Wick contrac¬ 
tions shown in Figures 23 c) and (d). 

The electromagnetic current operator in units of e has 
the flavor structure j® m = quU'y^u + qdd'y^d + qsS'y^s, 
and it is convenient to express this in terms of SU(3)f 
singlet and octet currents, 


f™= 7 §(?« + qd + qs) jl 

+ + q d - 2q s ) j*~ 

+ 72 (fe — Qd)jfj + i 

where 

il~ = (w7+ d'y^d + S 7 m s) 

jl~ = 75 («7 u.u + dj^d - 2 s 7m s) 

= 75 “ ^ d ) ’ 


(Dl) 


(D2) 


and where the ± superscript indicates the current’s 
G-parity. 

Considering Eq. |D2| in the SU (3 )p limit where the 
propagators for up, down and strange quarks are identi¬ 
cal (U = D = S = Q) 1 it is clear that the octet currents 
and cannot give rise to a disconnected contribu¬ 
tion of the type shown in Figure [23) I)), as the relevant 
part of the Wick contraction would be proportional to 


tr [ 7 fjU] - tr [ 7 fj.D] = tr [y^Q] - tr [ 7 M Q] = 0 
tr [y^U] + tr [ 7 m D] - 2 tr [y^S] = (1 + 1 - 2 )tr [y^Q\ = 0 , 
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respectively. In general, the singlet current could give rise 
to a disconnected contribution proportional to tr [y^Q], 
but since it enters the electromagnetic current with a 
weight of q u + qa + q s = +§ — | — § = 0 , it too does not 
contribute in practice. 

For the case of the form-factors of (I, I z ) = (1, +1) 
mesons like the 7 r + or p + , the initial and final state 
mesons are identical and have the same G-parity. The 
G-parity invariance of QCD (assuming isospin invari¬ 
ance) then ensures that the G-parity of the current must 
be positive and only can contribute. As such the ma¬ 


trix element is proportional to q u — qd = +§- 1 = +1 

as expected, and we see that there cannot be any discon¬ 
nected diagrams even away from the SU(3)p limil[^] 

For transitions like p + — > 7 r + 7 , where the G-parity 
flips, while there are disconnected contributions away 
from the SU(3)f limit, the logic presented above in¬ 
dicates that there are not in this calculation which is 
performed with exact SU(3)F symmetry. Only con¬ 
nected diagrams are required, and the matrix element 
is proportional to q u + qd- We note that in the SU(3)f 
case, the matrix element for K* + —> K + j is identical to 

p + —> 7T + 7. 
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